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ABSTRACT 


This  report  describes  a mathematical  model  for  the  ambient  noise  in  the  ocean 
caused  by  merchant  ships.  This  kind  of  noise  dominates  in  a frequency  range  in 
which  many  initial  detections  are  made  by  modern  sonar  equipment  and  hence  an 
understanding  of  it  is  of  considerable  interest.  ~ 

^A  stochastic  model  for  ambient  noise  is  defined  and  its  main  properties  are 
explored.  The  ambient  noise  process  is  a function  both  of  time  and  of  space, 
and  it  is  shown  to  be  a stationary  stochastic  process  as  a function  of  either  of 
these  variables.  Its  characteristic  function  is  calculated  as  are  its  first  two 
moments.  It  is  shown  that  the  ambient  noise  process  as  a function  of  time  has 
approximately  an  exponential  autocovariance  function  and  it  may  be  satisfactorily 
modeled  as  a Gauss-Markov  process.  Further,  it  is  shown  that  the  correlation 
between  the  ambient  noise  at  two  points  depends  to  a good  approximation  only  on 
the  propagation  loss  function  and  on  the  distance  between  these  two  points  (and  not 
materially  on  the  density  or  loudness  of  the  noise  field).  Finally,  an  effective 
and  practical  method  for  modeling  the  multidimensional  detection  process  is 
presented. 
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This  is  a report  by  Daniel  H.  Wagner,  Associates  to  the  Naval  Analysis 
Division  of  the  Office  of  Naval  Research  on  a research  investigation  under  ONR 
Contract  No.  N00014-76-C-0637.  The  report  is  addressed  to  the  modeling  of 
the  sonar  detection  process  when  there  are  multiple  sensors,  such  as  a field 
of  sonobuoys.  The  approach  taken  is  to  assume  that  the  random  component  of 
signal  excess  has  three  parts:  an  independent  component,  a correlated  component, 
and  a residual  component. 

It  is  hypothesized  that  the  correlated  component  is  principally  due  to 
ambient  noise,  and  a mathematical  model  for  ambient  noise  due  to  shipping  is 
defined  and  its  main  properties  discussed.  A convenient  method  for  simulating 
sample  paths  of  the  correlated  ambient  noise  fluctuations  is  then  given. 

We  would  like  to  acknowledge  the  cooperation  and  support  we  have  received 
during  this  work  from  Mr.  James  G.  Smith  and  Mr.  J.  Randolph  Simpson  of 
the  Naval  Analysis  Division  of  the  Office  of  Naval  Research. 
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SUMMARY 


In  sonar  detection  modeling  it  is  assumed  that  a detection  of  a target  occurs 
if 


x(t)  + i (t)  > 0 


at  some  time  t,  where  x(t)  = FOM(t)  - Nw(R(t))  is  the  deterministic  component  of 
signal  excess  at  the  sensor  at  time  t when  the  range  between  sensor  and  target  is 
R(t),  and  I (t)  is  its  random  component.  The  term  FOM(t)  represents  the  sonar 
figure-of-merit  at  time  t,  and  Nw  is  the  propagation  loss  function.  Various  models 
for  the  random  component  are  in  use;  these  and  the  available  results  in  computing 
cumulative  detection  probability  (cdp)  and  mean  time  to  detect  are  given  in 
reference  [a]. 

Now  suppose  there  are  n sensors,  and  n corresponding  signal  excess  histories 
xl*  • • • » ^n*  The  problem  is  to  hypothesize  a structure  tor  the  corresponding  ran- 
dom components  1 £n.  There  are  two  easy  ways  out: 

(1)  the  Ij's  are  all  independent; 

(2)  the  £ j's  are  completely  dependent,  i.  e. , § j = |2  = • • • = 4n- 

Another  possibility  is  to  compute  some  quantity  of  interest  such  as  cdp  or  mean  time 
to  detect  under  assumptions  (1)  and  (2)  and  average  the  results  (weighting  them 
equally,  as  in  reference  [b],  or  almost  equally,  as  in  reference  [c]).  Of  course, 
none  of  these  is  correct;  the  right  amount  of  correlation  is  surely  range  dependent, 
as  we  shall  see. 

We  will  hypothesize  that  each  £ { is  of  the  form 

4l(t)=4c(t)  + r,l(t)  + r1;(t) 

where  § is  identical  at  each  sensor,  Is  independent  of  r|£  for  all  i A and 
of  rik  for  allk,  r)  J is  independent  of  § c for  all  i,  rij  is  independent  of  £ c for  all  i, 


but  t| ^(t)  and  rik(s)  are  correlated.  For  these  reasons  £ may  be  called  the  residual 
component,  the  correlated  component,  and  tj!  the  independent  component.  The 
process  £c  represents  fluctuations  due  to  the  target,  rj  j represents  fluctuations  due 
to  ambient  noise,  and  rjj  represents  fluctuations  due  to  sensor  performance  and  propagation. 

In  Chapter  II  a mathematical  model  for  ambient  noise  due  to  merchant  shipping 
is  given.  It  is  assumed  that  the  ambient  noise  at  a point  P at  time  t is  given  by 


ri(P.t)  * 10  log  2j  Z f(R.(P,  t)) 

where  Z.  is  the  power  emitted  by  the  j^1  noise  source,  f is  the  power  loss  function 
and  Rj(P,t)  is  the  distance  between  the  j^1  noise  source  and  the  point  P at 
time  t.  The  identification  ^(t)  = ri  (P^  t)  - E ii  (Pj,  t),  where  Pf  is  the  position  of 
the  Ith  sensor,  will  be  made. 

The  characteristic  function  of  the  related  process  £ = lO^10  is  computed  as 
are  its  mean  and  variance.  It  is  shown  that  the  inclusion  of  the  effects  of  attenuation 
in  the  propagation  loss  function  may  be  crucial. 

In  Chapter  in  it  is  shown  that  rj  is  a stationary  stochastic  process  as  a function 
of  time.  It  is  also  argued  that  as  a function  of  t is  approximately  a Gauss-Markov 
process,  and  a parametric  estimate  for  the  relaxation  time  is  given. 

In  Chapter  IV  the  correlation  between  17  (P,  t)  and  t](Q,  t)  is  shown  to  be  ap- 
proximately independent  of  everything  but  the  propagation  loss  function  and  the  dis- 
tance between  P and  Q.  A method  for  estimating  ^(u),  the  correlation  in  ambient 
noise  between  two  points  separated  by  distance  u,  is  described  and  applied  to  a 
set  of  eight  environments.  The  qualitative  effects  of  convergence  zones  on  the 
function  are  described. 

Finally,  in  Chapter  V,  a method  for  simulating  the  correlated  processes  r^,  . . . , 
iln  is  given.  This  involves  simulating  n independent  copies  of  the  ti  process  and 
performing  a matrix  transformation  which  will  result  in  n processes  with  the  desired 
correlation  structure.  The  matrix  needed  for  this  transformation  turns  out  to  be 
the  square  root  of  the  matrix  (rik)))  > where  r ^ is  the  distance  between  sensor  i 
and  sensor  k. 

Our  model  for  ambient  noise  due  to  shipping  assumes  a homogeneous  density  of 
noise  sources.  That  is,  we  do  not  consider  cases  where  there  are  significant  de- 
partures from  homogeneity  such  as  heavily  traveled  shipping  lanes.  Most  of  our 
results  will  no  doubt  have  to  be  modified  in  order  to  treat  non -homogeneous  noise 
fields,  but  these  extensions  will  be  postponed  for  future  consideration. 
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THE  EFFECTS  OF  AMBIENT  NOISE  ON  THE  DETECTION 
PROCESS  AND  THE  MODELING  OF  DETECTION 
BY  A FIELD  OF  SENSORS 


CHAPTER  I 
INTRODUCTION 


In  a sonar  detection  model  one  hypothesizes  that  a sensor  is  in  contact  with 
I a target  at  time  t if 

i 

1 L 

x(t)  + £ (t)  > 0 


where  x is  a deterministic  function  (called  the  signal  excess)  obtained  from  the 
sonar  equation 


*(t)  - Ls<t)  - (LN(t>-NDI<t>)  - NR[|,t)  - Nw(R(t)) 
- FOM(t)  - Nw(R<t)) 


and  | is  a stationary  stochastic  process  representing  the  combined  uncertainties 
and  fluctuations  in  the  components  of  the  sonar  equation.  In  the  above,  FOM  is 
called  the  figure  of  merit,  Lg  is  radiated  noise,  LN  is  ambient  or  self  noise,  NDI 
is  the  directivity  index,  Njyj  is  recognition  differential,  R is  a function  giving 
distance  between  target  and  sensor,  and  Nw  is  the  propagation  loss  function.  The 
sonar  equation  and  its  components  are  discussed  fully  in  texts  such  as  reference  [d]. 


Given  a specific  mathematical  model  for  £ one  can  then  in  theory  compute  such 
items  of  interest  as  cumulative  detection  probability, 

cdp(t)  = Pr{x(s)  + £ (s)  > 0 for  some  s<  t}  , (1-1) 

mean  time  to  detect,  mean  time  in  contact,  etc.  See  references  [a]  or  [f]  for  surveys 
of  results  available  for  this  single-sensor  model,  and  reference  [e]  for  a discussion  of 
the  need  for  methods  of  estimating  cdp. 
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Id  this  report  we  discuss  a mathematical  model  for  that  component  of  £ which 
is  due  to  fluctuations  in  ambient  noise  due  to  shipping.  Thus  we  address  fluctuations 
contributed  by  the  function  Ljj,  and  then  only  when  L^j  is  dominated  by  con- 
tributions from  shipping  sources.  This  is  generally  in  the  50-500  Hz  frequency 
band,  an  interval  of  some  interest  in  sonar  detection.  Fluctuations  of  the  other 
components  may  be  modeled  by  the  methods  of  reference  [a]  until  such  time  as 
mathematical  models  for  them  are  available. 


This  ambient  noise  model  is  applied  to  derive  a method  for  treating  the  problem 
of  computing  cdp  for  a field  of  sensors.  That  is,  suppose  there  are  n sensors  and  a 

target's  behavior  produces  n signal  excess  histories  Xj and  n corresponding 

fluctuation  process  sample  paths  £ 1 £ n.  Then  we  are  often  interested  in 

computing  the  function 


cdp(t)  - Pr{x.(s)  + £ .(s)  > 0 for  some  s < t and  some  i = 1,  . . . , n}  , 


which  is  the  cumulative  probability  of  detection  by  some  member  of  the  field  of  n 
sensors.  (From  here  on  cdp  has  this  more  general  definition,  as  opposed  to 
(1-1). ) 

The  problem  in  computing  cdp  is  that  some  assumption  has  to  be  made  about 
the  joint  probability  distribution  of  the  £ j.  For  example,  if  £ i is  independent*  of  |k 
for  all  i^k,  then 


cdp(t)  = 1 - .Dj  (l  ~ cdp.(t)) , 


where 


cdp.(t)  = Pr{x.(s)  + £ .(s)  > 0 for  some  s < t} . 

On  the  other  hand  if  the  £ . are  completely  dependent  so  that  §=§„  = ...=§  , then 

i la  n 

cdp(t)  = Pr{  max  x.(s)  + £ . (s)  > 0 for  some  s < t} 

1 < i < n i l ~ 


The  process  is  independent  of  |k  if  the  random  variables  £j(s)  and  £k(t)  are 
independent  for  all  s and  t. 


which  is  of  the  form  (1-1)  with  x(s)  = maxj  < 4 < n x^s). 

In  reality  the  correlation  between  £^(t)  and  ^(t)  for  a given  t is  somewhere 
between  0 and  1.  That  is  because  each  sensor  is  subject  to  a variety  of  types  of 
fluctuation  some  of  which  are  shared  to  varying  degrees  by  the  other  sensors 
and  some  of  which  are  peculiar  to  the  sensor  itself.  We  hypothesize  that  each  £ ^ 
is  of  the  form 


ii<t)  = Sc(t)  + Tii<t)  + Tii<t) 


where  £ c is  identical  at  each  sensor,  r|  J is  independent  of  for  all  i and  of 
for  allk,  rj | is  independent  of  £ c for  all  i,  rjj  is  independent  of  £ c for  all  i,  but 
r|j(t)  and  r^s)  are  correlated.  For  these  reasons  £ c may  be  called  the  residual 
component,  q . the  correlated  component,  and  tj  I the  independent  component.  We 
envision  £ c as  deriving  from  fluctuations  due  to  the  target,  ri ! as  deriving  primarily 
from  fluctuations  in  Njjj,  N^,  and  propagation  loss  anomalies  in  the  immediate 
vicinity  of  the  sensor,  and  as  deriving  primarily  from  fluctuations  in  the  ambient 
noise  field.  As  such  -rj  j and  qjj  will  be  positively  correlated  since  they  share  the 
same  ambient  noise  field,  albeit  a different  view  of  it. 

Now  suppose  through  theory  or  experiment  one  has  an  estimate  of  p^,  the  cor- 
relation between  rj . (t)  and  r^(t)  for  each  t > 0.  How  does  one  go  from  the  matrix 
((  pik))  correlations  to  a method  for  simulating  sample  paths  r)  j , . . . , r|  ? This 

question  is  resolved  in  the  final  chapter.  In  the  intervening  chapters  we  define  the 
ambient  noise  process,  investigate  its  elementary  properties,  show  that  it  is  a 
stationary  stochastic  process  and  that  it  can  be  adequately  approximated  as  a 
Gauss-Markov  process. 

It  should  be  mentioned  here  that  methods  for  computing  cdp  analytically  are  not 
available  except  for  the  cases  where  the  £ j are  all  independent  or  completely 
dependent.  For  the  more  realistic  case  of  correlated  £ j's,  Monte  Carlo  methods 
must  be  used. 


THE  AMBIENT  NOISE  MODEL 


It  ceased;  yet  still  the  sails  made  on 
A pleasant  noise  till  noon, 

A noise  like  of  a hidden  brook 
In  the  leafy  month  of  June, 

That  to  the  sleeping  woods  all  night 
Singeth  a quiet  tune. 


Samuel  Taylor  Coleridge, 
The  Rime  of  the  Ancient  Mariner 


Suppose  that  at  time  t,  ambient  noise  sources  to  be  thought  of  as  ships  are 
located  at  points  { Sj  (t)}  j in  the  Euclidean  plane  where  J is  a countable  index 
set.  The  ship  is  assumed  to  be  an  omnidirectional  noise  source  producing 
a noise  level  of  Nj  db  (referenced  to  some  fixed  intensity,  at  some  fixed  frequency 
of  interest).  The  ship  is  assumed  to  be  on  a fixed  course  0^  at  a constant 
speed  V: . 


Then  for  each  point  in  the  plane  and  each  time  t there  is  defined  a random 


variable 


where  d(P,  Q)  is  the  distance  between  the  points  P and  Q.  (In  all  that  follows 
log  will  denote  common  logarithms  and  In  will  denote  natural  logarithms. ) The 
quantity  r|(P,  t)  is  the  magnitude,  in  db,  of  noise  reaching  the  point  P at  time  t, 
and  t)  is  a stochastic  process  as  a function  of  P and  t.  It  will  sometimes  be  referred 
to  as  t|  (P)  when  t is  fixed,  or  as  ri  (t)  when  P is  fixed.  To  relate  it  to  the  previous 
chapter,  rj - (t)  will  be  defined  as  ^(Pj,  t) -E  ri(Pj, t)  where  Pj  is  the  location  of  the  i*“ 
sensor. 


We  will  also  have  need  of  the  power  process  f defined  by 


£(p.t) 


sJio(Nj-Nw(d<p-sj<t»))  /10 


,10T,(P,t)/10_ 


The  initial  positions  Sj(0)  of  the  noise  sources  at  time  t = 0 will  be  assumed 
chosen  in  one  of  two  ways.  For  the  Poisson  noise  process  we  assume  that 
(sj(°)}  j€  j ifi  a sample  from  a 2-dimensional  Poisson  process  with  intensity  \ > 0. 
For  the  bounded  noise  process  we  assume  that  some  fixed  finite  number  M > 0 of 
points  {8,(0)}.“  are  selected  independently  and  according  to  the  uniform  distribution 
on  some  bounded  measurable  subset^ of  the  plane. 

The  initial  courses  6,  will  be  assumed  drawn  independently  from  the  uniform 
distribution  on  [0,  2tt).  These  courses  remain  constant  for  a Poisson  noise  process. 
For  a bounded  noise  process  the  noise  sources  are  assumed  to  reflect  off  the 
boundaries  of  when  they  are  run  into.  The  speeds  V:  and  noise  levels 
Nj  are  assumed  drawn  independently  from  distributions  Fy  and  FN,  respectively. 
Both  vj  and  Nj  remain  constant  for  all  t > 0. 

I 

Before  proceeding  to  describe  the  probability  distribution  of  p (P,  t)  for  fixed 
P and  t,  we  need  to  mention  some  properties  of  the  spatial  Poisson  process. 

Suppose  { Sj}  je  j is  a sample  from  a spatial  Poisso)  process  with  parameter  X, 
the  mean  number  of  points  per  unit  area.  If  A is  a measurable  subset  of  the  plane 
let  n (A)  denote  its  area  and  let  \ (A)  be  the  number  of  points  from  the  set  { S:}  j 
that  are  in  A.  Then  the  following  properties  will  be  useful: 

(i)  The  random  variable  x (A)  has  a Poisson  distribution  with 
parameter  (A),  that  is, 

Pr{  X (A)  = k}  = (A)  -^:(A))  , k = 0,  1 

(ii)  If  Aj  and  A„  are  disjoint,  then  x(Aj)  and  x(^2^  are  independent 
random  variables. 

(iii)  Given  that  x(A)  = k,  then  the  k points  in  A are  distributed 
independently  and  according  to  a uniform  distribution  on  A. 

To  derive  the  characteristic  function  <p  of  £(P,t)  we  note  that  P can  be  taken 
as  the  origin,  without  loss  of  generality,  and  that  f is  of  the  form 


‘-j&W 

where 
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■ m 


r T f 

1 i 

L 


Zj  - ioY10, 


r is  the  distance  between  S.(t)  and  the  origin,  and 
j J 

f,r)  - X0-Nw'r»/10 


is  called  the  power  loss  factor  at  range  r. 

Theorem.  Let  { Zj}  j€  j be  a sequence  of  independent  non-negative  random 
variables  with  common  distribution  function  Fz,  let  {S:}  j€  j be  a sample  from  a 
spatial  Poisson  process  with  parameter  \ and  let  r.  be  the  distance  of  S.  from  the 
origin,  for  je  J.  Then  J ^ 


£ = .2  Z.  f(r.) 
je  J J ' J 


has  the  characteristic  function  cp  given  by 


<p(y)  = exp{2irX/“  /“  (elzyf(r)  - 1)  r dr  d Fz(z)}  , 


Proof*  Let  £m  denote  the  random  variable 


£ = 2 Z.  f(r.) 

m rj  < m J J 

with  characteristic  function  <pm>  for  m = l,  2,  ...  . Then  {£  } converges  a.  s. 
to  £ so  that  by  the  dominated  convergence  theorem  cp  (y)  converges  to  (p  (y)  for 
all  y.  Let  Bm  denote  the  disc  with  radius  m,  centered  at  the  origin.  Then 

( p (y)  = Eeiy^m 
mw 

= Jo  E(eiyf  m I X (Bm)  = k)  Pr{  x (Bm)  = k} 

-kl0  E(elyfn,|  x(Bm)  . e-^m? 

I 

using  the  fact  that  £ , conditioned  on  X(Bm)=k>  is  the  sum  of  k independent  and 
identically  distributed  terms.  Thus 
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*«w  - A C ‘!,yZ'<^,  «*,»!  k e-X ’m2 

K- 

-e-^“2  exp{2„*/;/0m  e^«r»rdrdFz(z)} 

-exp{2,A;;/”(el>’zW-l)rdrdFz,z,} 

- exp{27rX  /“  /”  (eiyzf(r)  - l)  r dr  d Fz(Z)} 

a8®-®,  and  the  theorem  is  proved.  In  the  following  assume  that  the  Z-  are  bounded. 
Corollary  1.  The  random  variable  £ is  finite  a.  s.  if  and  only  if 


fQ  r f (r)  dr  < ® . 


(II- 2 ) 


Proof-  From  (II-l)  it  follows  that  £ < ® a.  s.  if  and  only  if 


/“  (eiyzf(r)  - l)rdr  < <°. 


(II- 3) 


iyzf(r) 


l)  r = iyzf(r)  r + r - o^f(r)^  , r - ® , 


so  that  (U-3)  holds  if  and  only  if  (II— 2)  holds. 

Corollary  2.  If  £ < « a.  s.  and  p z = j“  z d Fz(z  ),then 

E £ =2*xHzfo  rf<r)dr  < ®» 


(II- 4) 


= 2t rX  E(Z2) / ' r f(r)2  dr  < 


Proof-  This  follows  from  (II-l)  and  the  well-known  formula  E £k  =i"k^<k)( 
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Example.  Suppose  Nw  is  of  the  form 


Nw(r)  = 1 0(a  log  r + /3r),  r > rQ  (in  yds), 


so  that 


f(r)  = r~a  lO~Pr,  r > rQ. 


(II- 5) 


The  parameter  at  is  called  the  spreading  constant  and  10 p is  che  attenuation 
coefficient.  Then  Corollary  1 implies  that  ambient  noise  is  finite  if  and  only  if 
/3  > 0 or  a > 2.  In  particular,  if  one  is  modeling  ambient  noise  and  assuming 
spherical  spreading  or  less  (a  < 2)  then  one  must  consider  the  effects  of  attenuation 
(i.  e. , assume  /3  /0)  or  else  spurious  results  will  be  obtained.  This  effect 
might  not  have  been  anticipated,  especially  for  low  frequency  propagation,  since 
attenuation  at  low  frequencies  is  apparently  insignificant.  For  example,  at- 
tenuation is  only  about . 0005  dbAyd  at  50  Hz,  and  yet  must  be  included  in  one's 
model  if  one  is  assuming  an  idealized  infinite  planar  ocean.  Of  course  if  one  is 
using  a bounded  noise  process  and  the  region  has  diameter  say  3,  000  miles,  then 
attenuation  can  be  safely  ignored. 

The  mean  value  given  by  (H-4)  can  be  evaluated  in  closed  form  for  some  cases 
of  (II-5).  Suppose  f is  of  the  form  (II-5)  with  rQ  = l yd.  Then 

,°°  ,,  . , _ , , -<D  -a  +1  -lnl0/3r, 

fQ  f(r) r dr  = f(r)rdr  =/  r e H dr 


= (lnl0/3)"'2r(-o+2,  lnl0/3), 


where  T is  the  incomplete  gamma  function,  and  the  latter  can  be  evaluated  ex- 
plicitly for  o=l,  3/2,  2.  Table  II— 1 gives  some  sample  calculations. 

At  this  point  we  should  mention  something  about  the  need  for  having  two 
mathematical  models:  the  Poisson  noise  process  and  the  bounded  noise  process. 
The  main  reason  is  that  the  Poisson  noise  process  is  much  more  suitable  for 
proving  theoretical  results,  but  we  need  the  bounded  process  to  represent  the  real 
world.  Moreover  to  simulate  a sample  path  of  the  time  series  t-*ri(P,  t)  we  can 
only  generate  a finite  number  of  noise  sources  so  the  bounded  noise  process  is 
required.  Many  of  our  conclusions  are  based  on  an  examination  of  one  of  these 
processes  and  we  will  often  glide  effortlessly  to  the  corresponding  conclusion  for 
the  other  process,  when  it  seems  plausible.  t 


TABLE  II-l 


TYPICAL  VARIATIONS  IN  AMBIENT  NOISE 
AS  A FUNCTION  OF  SPREADING  CONSTANT 


THE  AMBIENT  NOISE  TIME  SERIES 


In  this  chapter  we  etudy  properties  of  the  time  series  t-*Ti(P,  t),  for  a fixed  P. 
Thus  we  suppress  the  dependence  on  P and  write 


nW  = 10  log  ^ io^J-VYDV10 

= 10  logZ.Z.  f(r.(t)), 


where  r-  (t)  is  the  distance  of  Sj (t)  from  some  fixed  point  P.  In  Chapter  IV  we 
will  study  the  spatial  function  P— r](P,  t),  with  t fixed.  Recall  that  Sj(t)  is  the 
position  of  the  j**1  noise  source  at  time  t.  If  t]  is  a Poisson  noise  process,  then 
the  choice  of  P is  immaterial  and  we  assume  it  is  the  origin.  If  ^ is  a bounded 
noise  process  then  the  choice  of  Pe^  is  relevant  to  the  statistical  behavior  of  the 
process  t-*ri(P,t). 

We  will  show  that  a bounded  noise  process  on  a rectangle  or  a Poisson 
noise  process  is  a stationary  stochastic  process.  It  is  conjectured  that  an  arbitrary 
bounded  noise  process  is  a stationary  process  but  we  have  not  been  able  to  prove 
this.  Also  we  discuss  the  sample  path  behavior  of  t-*-r|(t),  showing  that  this  time 
series  is  approximately  a Gauss-Markov  process  and  giving  a formula  to  approximate 
the  relaxation  time. 

Theorem.  The  bounded  noise  process  t -*  rj  (t)  on  a rectangle  is  a stationary 
stochastic  process. 

Proof.  We  must  show  that  the  distribution  of  the  random  vector  (t)  (tj+A),  . . . , 

r|  (tj^+A))  is  independent  of  A for  any  integer  k,  and  time  points  tj t^.  Since 

r|  (t)  depends  deterministically  on  the  random  variables  Sj(t) S^(t),  Nj,  .... 

Nm,  clearly  it  is  sufficient  to  show  that  the  distribution  of  the  kM -dimensional  vector 


^®j(tj+A),  ...»  Sj(t^+A),  . 


SM(tl+A)’  •••'  SM(tk+A)) 


*-  ! 


fc:  « 


is  independent  of  A.  Moreover,  since  Sj(tj+A)  is  independent  of  S.(t  +A)  for 
any  i^j,  it  suffices  to  show  that  the  distribution  of 


(s^+A) S^+A)) 


is  independent  of  A. 


First  we  show  that  the  distribution  of  S^t)  is  uniform  on  for  any  t > 0. 
Suppose  without  loss  of  generality  that  the  vertices  of  ^5  are  at  (0, 0),  (a,  0),  (a,  b), 
(0, b),  a>b.  Let  h:(x,y)-h(x,y)  be  the  probability  density  of  Sjft),  (x,y)e^. 

Let  h(* , • I Vj)  denote  the  conditional  density  of  Sj(t),  conditioned  on  the  velocity 
Vj.  Then 


h(x,  y)  = J0  h(x,  y | Vj)  d Fv(Vl), 


so  it  suffices  to  prove  that  h(x, y|  v^)  is  uniform  on  for  a given  fixed  v..  Suppose 
at  first  that  t < b/(2  v^).  We  have 


h(x,y|v1)  = h(x,y|vr  flj)  dflj 


where  6 is  the  initial  heading  of  this  noise  source.  Thus 

h(».y| vi)  - -27 1(;0'/2  ♦ ;;/2 ♦ /2'/2 + 4*/2) k<*.y I v2. 


= W'Jo^2  h(x’ylvr  ei)  dei  + Jq  /2  h(x,yl  V 7r“e1)d01 


+ r0/2  h(x»y|  Vj,  flj+Tr)  dex  +f* 72  h(x, y | Vj,  dej 


tt/2 


= “/q /2  T{h(X,ylVl’  V + h(x,y|v1,  ir-flj) 


h(x,y|  Vj,  Cj  + jt) +h(x,y|v1,  -©j)!  dfl^ 


A glance  at  Figure  III-l  illustrates  that  the  integrand  in  the  last  expression  above 
is  equal  to  1/ab  for  any  t < b/(2Vj). 


) 

L 


-12- 


Note:  Number  in  each  box  is  the  value  of  the 
resulting  density  xab. 


v t cos  6 


Thus  S| (t),  conditioned  on  the  value  of  vj,  is  distributed  uniformly  for  t suf- 
ficiently small  and  hence  for  any  t,  since  the  effect  of  an  arbitrary  t can  be  broken 
up  into  successive  sufficiently  small  steps.  Thus  Sj(t)  is  distributed  uniformly, 

independent  of  t.  Finally  the  joint  distribution  of  (S^t.+A),  S.^+A) Sj(t+A)) 

depends  only  on  whether  the  given  points  can  be  joined  oy  an  admissible  path  and  an 
the  speed  Vj  required  to  traverse  this  path,  and  not  on  A.  This  completes  the  proof. 


The  details  of  the  proof  make  clear  that  to  extend  this  result  to  an  arbitrary 
bounded  noise  process  one  need  only  show  that  our  model  for  motion  preserves  the 
uniform  distribution  on  the  region^.  Our  method  of  proving  this  for  a rectangular 
^obviously  will  not  apply  to  a general  region  which  does  not  have  the  symmetries 
of  a rectangle. 

The  fact  that  a Poisson  noise  process  is  stationary  is  obvious.  The  points 
{Sj(t)}je  j are  a sample  from  a spatial  Poisson  process  with  intensity  \ for  any 
t > 0,  and  the  joint  distribution  of  the  (infinite)  random  vectors 

fsjVA,t)£J' {s)<VA,1)cj 

depends  only  on  whether  Sjftj+A),  ....  Sj(t +A)  can  be  joined  by  a path  consistent 
with  our  model  for  motion,  for  all  je  J.  This  in  turn  depends  on  the  quantities 
(tj  +A)  - (tj  _2+A)  = tj  -tj  _j  but  not  on  A. 


We  now  turn  the  discussion  to  the  sample  path  behavior  of  an  ambient  noise 
process.  It  turns  out  that  the  sample  paths  of  a bounded  noise  process  have  as 
many  derivatives  as  the  propagation  loss  function  Nw,  except  that  a bounded  noise 
process  will  have  isolated  non-differentiabilities  corresponding  to  time  points  at 
which  noise  sources  bounce  off  the  boundaries.  This  is  proved  by  noting  that  the 
sample  path  is  the  logarithm  of  the  sum  of  a finite  number  of  functions  each  of  which 
has  the  required  smoothness.  If  Nw  is  differentiable,  / rf'(r)dr  < ®,  and  both 
Zj  and  Vj  have  finite  means,  then  the  corresponding  Poisson  noise  process  sample 
paths  are  differentiable  with  no  exceptions  since  the  noise  sources  do  not  bounce 
off  the  boundaries.  This  follows  from  examining  the  difference  quotient 

h"1  (£  (t+h)  -£(t))  =h_1  2 Z.jj(r.(t+h))  - f (r.(t))J 
- 2 Z.  f'(ryt))  r!(t)  as  h - 0. 

Figure  in-2  shows  a sample  path  for  a bounded  noise  process  on  a square 
region  with  area  A = 9x10®  square  miles,  with  M = 900  noise  sources,  the  v's 
uniformly  distributed  between  15  and  25  knots,  N;  = 180  db//l^  Pa//1  Hz,  and 
propagation  loss  1A  from  Chapter  IV.  The  four  large  excursions  during  this  eight 
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£2 


Time  (days) 


day  period  correspond  to  time  periods  when  a noise  source  comes  close  to  the 
sensor.  In  each  of  these  four  cases  a noise  source  comes  within  6 nm  of  the 
sensor. * 

Since  h is  a stationary  stochastic  process  it  has  an  autocorrelation  function 
* and  it  is  of  great  interest  to  compute  *.  Unfortunately,  analytic  methods  have 
so  far  failed  to  yield  any  results  even  for  the  simplest  Nw>  or  for  the  simpler 
power  process  £ . Some  Monte  Carlo  estimates  have  been  performed  for  the  case 
which  yielded  the  sample  time  series  in  Figure  III-2,  with  variations  on  M,  how- 
ever, M = 900,  300,  100,  and  40.  ** 

Figure  III— 3 shows  the  estimated  autocorrelation  functions  plotted  on  semi- 
log paper.  We  note  that  the  data  points  indicate  an  almost  linear  relationship, 
i.  e. , ^(t)  is  approximately  of  the  form  e_<5t.  Linear  fits  for  the  parameter  6 
have  been  made  and  tabulated  in  Table  IH-1. 


TABLE  III-l 

ESTIMATED  RELAXATION  TIMES 


Number  of  ambient 
noise  sources,  M 

Estimated  Relaxation 
Time  1/6  (hrs. ) 

v 1 VA/M 

Ratio 

900 

3.61 

5.00 

.72 

300 

5.41 

8.  70 

.62 

100 

8.47 

15. 15 

. 56 

40 

15.38 

23.  80 

.65 

At  any  given  time  one  of  the  noise  sources  is  closest  to  the  fixed  point  P.  Con- 
sider the  mean  time  v that  such  a source  persists  as  the  nearest  neighbor  to  P. 

It  might  be  conjectured  that  the  "relaxation  time"  of  the  process  (the  coefficient 


* In  most  naval  applications  such  near  field  noise  sources  are  not  considered 
part  of  the  ambient  noise  since  they  are  "identifiable.  " There  is  a problem 
with  this  viewpoint,  however,  since  so  many  of  the  noise  sources  are  potentially 
identifiable  if  only  the  effort  were  to  be  exerted.  In  any  event  a definition  such 
as  ours  is  necessary  in  order  to  have  a mathematically  tractable  model. 

**  For  a fixed  one  need  only  vary  M (or  X for  a Poisson  noise  process)  to  see 
all  the  possibilities.  A change  in  the  mean  level  of  the  source  levels  N*  cor- 
responds to  a shift  in  the  mean  level  of  t).  A multiplicative  change  to  the  dis- 
tribution of  the  Vj  corresponds  to  a multiplicative  change  in  the  time  scale  of  tj. 
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Time  (hrs) 


M - 100  U_ 


FIGURE  III -3 

ESTIMATED  AUTOCORRELATION  FUNCTIONS  FOR  THE  AMBIENT  NOISE 
PROCESS: VARIOUS  SHIP  DENSITIES 


M 900 


M = 300 


1/6  In  the  autocorrelation  function  exp(-6  |t|))  is  somehow  related  to  v.  This  is 
because  the  nearest  point  to  P exerts  the  strongest  influence  on  q and  because, 
on  the  average,  q comes  under  a new  influence  every  v units  of  time.  If  the  M 
noise  sources  were  spread  out  evenly  across  the  region^?  and  the  particles  all 
moved  in  the  same  direction  at  the  mean  velocity  V = 20  knots,  then  nearest 
neighbor  changes  would  occur  every  v'VA/M'  hours.  This  quantity  has  been 
calculated  for  v = 20  knots,  /X  = 3000nm,  and  M = 900,  300,  100,  and  40,  and  entered 
in  Table  in-1  also.  We  notice  the  rather  good  agreement  between  estimated 
relaxation  time  and  the  quantity  (2/3)  7 * VA/M.  For  a Poisson  noise  process  the 
corresponding  quantity  to  use  would  be  (2/3)  v-1  A.-*/2. 

It  is  known  (Doob's  Theorem,  reference  [g],  p.  321)  that  a stationary  Gaussian 
process  is  a Markov  process  if  and  only  if  its  autocorrelation  function  is  of  the 
form  exp(-6  |t| ).  Now  in  our  case  we  have  a process  q which  is  stationary  (shown 
in  this  chapter),  which  is  approximately  Gaussian  (shown  by  using  a central  limit 
theorem  for  power  sums  in  reference  (hj),  and  which  has  an  autocorrelation  function 
which  is  approximately  of  the  form  exp(-6  |t| ) (shown  in  this  chapter).  Therefore  we 
are  encouraged  by  Doob's  Theorem  to  conclude  that  q is  approximately  Markovian 
and  hence,  using  its  other  properties,  that  q is  approximately  a Gauss-Markov 
process. 

Thus  we  have  concluded  that  for  any  fixed  P the  time  series  t— q(P,  t)  is  ap- 
proximately a Gauss-Markov  process  with  autocorrelation  function 

*(t)  = exp(-  J-  vVT  tj  (III-l) 

or 

*(t)  = exp^-  |-  iJf  t^.  (HI-2) 

This  conclusion  is  of  extreme  importance.  It  allows  us  to  simulate  the  time 
series  t-*q(P,  t)  (and  hence  the  qj  in  which  we  are  ultimately  interested)  without 
simulating  the  field  of  noise  sources  { Sj}  j£  j which  would  be  very  cumbersome 
and  time  consuming.  We  will  return  to  this  point  in  Chapter  V. 


1 


CHAPTER  IV 


SPATIAL  CORRELATION  OF  AMBIENT  NOISE 


In  this  chapter  we  look  at  the  process  P-*ti(P,  t)  at  a fixed  time  t, 
,1<P>=10  1ogSj10^Vrj<P>>)/10 


= 10  log  2jzj  f(r.(P» 


= 10  log  £<P), 


where  Tj(P)  is  the  distance  between  the  noise  source  and  the  point  P,  at  some 
fixed  time  t. 


For  a Poisson  noise  process  tj  is  a stationary  process  as  a function  of  P,  but 
no  bounded  noise  process  can  be  stationary.  For  the  rest  of  this  chapter  we 
assume  tj  is  a Poisson  noise  process. 


I 4 


: vr 


An  analytic  expression  for  the  correlation  between  £ (P)  and  £(Q)  can  be  given. 

/ 

Theorem.  Let  P = (0,  0)  and  Q = (u,  0).  Then  the  correlation  p between  £ (P) 
and  £ (Q)  is  given  by 

*2tt  /*°° 


rj  r f(r)  f ((r2  + u2  - 2ru  cos  & dr  d 6 

P ~ »ao 


2t r f 

J 0 


r f(r)  dr 


(IV-1) 


Proof.  Let  rj  and  rj  be  the  distances  of  the  j**1  point  (x.,  y ) of  the  Poisson 
pie  from  the  points  P and  Q respectively.  Let  Bn  be  the  disk  with  radius 


sample 

n centered  at  the  origin.  Let 


► ■*< , 


£ = E Z,  f(r  ) 
(Xj,yj)eBn  J J 
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n 


(Xj,yj)€Bn  J i 


so  that  £n~*  £ (P)  a.  s.  and  — £ (Q)  a.  s. , and  hence 


E fn  £n  "■  E £(P)  £(Qh  as  n -*00- 


Lex  xn  be  the  number  of  points  from  the  Poisson  sample  in  the  set  Bn. 
Xn  is  a random  variable  with  mean  A7rn2.  Let  JR  be  the  subset  of  indices  j 
that  (Xj , y j ) e BQ.  Then 


= E E Z-  f(r.)  f(r!) 
jeJn  } } ) 


+ EE  zi  E Z ; f(rl ) 

jeJn  J J UJn 

l 

•27&*n 


2 I . o 2 lv 

= E(x  )E(Z)  -/  I f(r)  f((r  +u  -2ru  cos  6)2)  r dr  dfl 

n 7rn“  •'0  */0  v ' 


CO 

+k?lE 


E z f(r  ) £ -Z,  f(r}) 
jeJ  3 J /Tj  ‘ 1 


leJr 

l /j 


X = k 
*n 


Pr(xn  = k) 


-A  E(Z2)^jT  f(r)  f((r2  +u2  -2r  cos  0)^)  r dr  d 0 

+ k=l  k(k_l)  Pr(xn=k)  (EZ)2^-E^  f(r)  r dr^ 
/'27r/“n  i 

x ~2  Jq  Jq  f((r2  +u2 -2ru  cos  0)2)  r dr  d0 
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Then 

such 


Here  we  have  used  the  fact  that  given  there  are  x points  in  B , they  are  distributed 
independently  and  uniformly.  The  formula  (IV-1)  now  follows  since  the  bracketed 
quantity  above  is  none  other  than  E f (P)  and  since  a2  p = 2tt\  E(Z2)/”  r f(r)2  dr. 


We  note  that  correlation  is  independent  of  X and  of  the  Z distribution,  depending 
only  on  the  propagation  loss. 

An  approximate  formula  for  the  correlation  ^(u)  between  t|(P)  and  r|(Q),  P and  Q 
separated  by  distance  u,  can  be  given  which  we  will  use  to  show  that  the  function 
is  also  approximately  independent  of  X and  the  Z distribution. 


Suppose  X and  Y have  a bivariate  Gaussian  distribution,  each  with  mean  m and 
variance  t2  and  correlation  p q.  Then  elementary  computations  show  that  the 
variables  io^/IO  and  10Y/10  have  mean  n,  variance  a2,  and  correlation  p given  by 


(IV -2) 


_ ln(l  + pg2/p2) 

° ln(l+a2/M2) 

But,  as  has  already  been  noted,  tj  is  approximately  Gaussian  and  so  the 
correlation  between  t]  (P)  and  t](Q)  is  given  by  (IV -2)  where  p is  the  correlation 
between  £ (P)  and  £ (Q)  and  is  given  in  (IV-1).  The  quantity  a2/p2  appearing  in 
(IV -2)  is  given  by 

_ E(Z2)  £ t<r)2  r dr 

9 9 oo  9 ' v1  v 

p 27rX(EZ)  (f ^ f(r)  r dr) 

using  Corollary  2 of  Chapter  I.  If  f is  of  the  form  (II— 5)  with  a = 1,  for  example*, 
then 


„ E(Z2) 

2ir  A (EZ)2 

= x-i  E(Z2) 
(EZ)2 


(lnlO)2  /32  In  (2/J) 
x 5 x 10-6, 


(IV -4) 


for  (3  = . 001  db/nm.  Now  for  reasonable  values  of  A and  E(Z2)/(EZ)2  such  as 
A > . 0001  and  E(Z2)/(EZ)2  < 2,  we  will  have  o2/p2  <<  1,  and 


Thus,  as  asserted,  the  correlation  between  t|(P)  and  ti(Q)  depends,  to  a good 
approximation,  only  on  propagation  loss. 

This  fact  is  extremely  useful,  since  it  allows  us  to  make  predictions  concerning 
spatial  correlation  of  ambient  noise  based  only  on  propagation  loss  in  an  area,  with 
no  need  for  predictions  on  the  density  or  noisiness  of  the  traffic.  If  a Monte  Carlo 
simulation  is  employed  one  need  only  plug  in  some  density  A and  constant  noise  Z 
and  estimate  correlation  for  a given  f;  there  is  no  need  to  vary  A and  to  see  the 
effects  of  such  variation. 


The  only  value  of  a for  which  the  calculation  (IV-3)  can  be  performed  in  closed 
form  is  o =1,  and  hopefully  the  approximations  to  follow  are  still  close  for 
1 <a  <2. 


■BBMDH 
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In  principle  (IV-1)  can  be  computed  by  numerical  analysis  and  the  approximation 
(IV-2)  employed,  with  a2/p2  given  by  (IV-3).  However,  it  is  easier  to  estimate 
as  a function  of  spatial  separation  via  a simulation,  and  thereby  relieve  any  mis- 
givings about  the  approximation  (IV-2),  and  avoid  the  difficult  numerical  analysis 
involved  in  (IV-1). 

Figures  IV-1  to  IV -4  show  some  spatial  correlation  functions  ^ corresponding 
to  the  potpourri  of  propagation  loss  functions  graphed  in  Figures  IV-5  to  IV-12. 

In  each  case  a density  of  . 0001  ships/(nm)2  with  source  level  180  db  was  assumed. 
For  ranges  past  250  nm  propagation  loss  was  assumed  to  be  of  the  form 
b + 12. 21og(2000r),  r in  nm,  the  constant  b chosen  to  make  propagation  loss  con- 
tinuous at  r = 250  nm.  In  each  case  ^(u)  was  estimated  for  14  values  of  u ranging 
from  1 to  75  nm,  1000  replications  in  each  case. 

The  main  conclusions  noted  from  an  examination  of  these  examples  are: 

(1)  correlation  is  significant  (say  greater  than  . 3)  out  to  a range, 
depending  on  propagation  loss,  that  may  be  anywhere  from 
10  to  100  nm, 

(2)  convergence  zones  induce  ripples  in  the  correlation  function, 
and, 

(3)  CZ’s  induce  a psuedo-randomness  in  the  sense  that  if  the 
oscillations  in  propagation  loss  caused  by  CZ's  have  large 
amplitude,  then  correlation  goes  to  zero  quickly. 

The  last  point  makes  sense  upon  reflection:  a slight  change  in  the  position  of 
a noise  source  may  alter  its  effects  on  a sensor  from  something  significant  to 
something  insignificant  and  vice  versa.  Thus  the  amount  of  coherence  in  the  non  e 
experienced  at  two  points  becomes  less  predictable,  i.e. , correlation  goes  down. 


CORRELATION  FUNCTIONS  CORRESPONDING  TO 
ENVIRONMENTS  1A  AND  IB 
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CHAPTER  V 


MODELING  MULTIPLE  SENSORS 

Finally  we  turn  to  the  central  problem  of  this  report,  that  is,  modeling  the 
detection  process  for  a field  of  sensors.  We  suppose  there  are  n sensors  and 
that  all  questions  concerning  detection  are  latent  in  the  functions  Xj  + lj  where  x. 
is  the  deterministic  component  of  signal  excess  and  £ . is  its  stochastic  component 
with  1 

I j(t)  = £c(t)  + Tij(t)  + Tij(t),  for  t > 0,  i = 1 n. 

We  assume  that  £c  is  independent  of  each  and  each  q.\  t).  is  independent  of 

% for  all  i andk,  and  tj'  is  independent  of  for  all  i/k.  T^ie  components  £ , 

ti  j,  ....  may  be  modeled  using  the  methods  of  reference  [a],  that  is,  as  (X,  a)- 
jump  processes  or  Gauss-Markov  processes  or  combinations  thereof. 

In  previous  chapters  we  have  introduced  a model  for  the  correlated  components 

’ll This  is  done  by  defining  a spatial  and  temporal  process  ti  and  taking 

i (t)  = r|  (Pj,  t)  - E r|  (Pj,  t)  where  Pj  is  the  location  of  the  i^  sensor.  We  showed, 
moreover,  that  each  is  approximately  a Gauss-Markov  process  with  autocorrelatioi 
function  ^ given  by  (III— 1)  or  (in-2).  This  is  extremely  fortuituous  because  the  ti 
process  is  very  cumbersome  to  simulate,  while  the  Gauss-Markov  process  can  be 
simulated  using  a simple  recursive  procedure  which  will  be  described  below. 

Our  object  now  is  to  give  a procedure  for  simulating  the  processes  r^,  ....  ^ 
such  that  each  process  t)  is  a Gauss-Markov  process  with  the  required  autocorrelation 
function  * and  such  that  ^.(t)  and  r)^(t)  have  correlation 

Pik  =<^(rik) 

where  r^  is  the  distance  between  sensor  i and  k and  ^ is  the  spatial  correlation 
function  discussed  in  Chapter  IV.  (It  should  be  noted  that  we  do  not  have  to 
specify  the  correlation  between  ti  . (t)  and  qk(s)  for  t /s,  since  for  a Gauss-Markov 
process  this  correlation  is  already  determined  by  specifying  pik  and  *,  and  is  equal  j 

to  *<t-s).  This  is  a trivial  consequence  of  the  recursive  relationship  (V-l)  below. ) 
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Before  proceeding  it  is  necessary  to  mention  some  additional  concepts. 
An  nxn  matrix  of  real  numbers  A = ((a^))  is  said  to  be  nonnegative  definite  if 

A A ”u  *]  >i  i ° 


for  all  real  numbers  y . 

1 n 

If  X-  (X- Xjj)  is  an  n-dimensional  random  vector  each  of  whose  com- 

ponents has  finite  variance  then  the  nxn  matrix  whose  (i,  j)  component  is  given  by 


E(Xt  - EX^fX.  - EX.) 

is  called  the  covariance  matrix  of  X.  The  matrix  with  (i,  j)  component 


E(Xj  - EX.)(X.  - EX.)/a.o. 


is  called  the  correlation  matrix  of  X,  where  a.  is  the  standard  deviation  of  X,.. 
A correlation  matrix  is  also  a covariance  matrix.  V 

A covariance  matrix  A = ((ai,))  of  a random  vector  X is  nonnegative  definite 
because  for  any  yp  . . . , yn 


n n 

A A \i  yj  yi 

- e(A  yi<v  exi>)2  ^ »• 


If  A is  a nonnegative  definite  matrix  then  there  exists  a lower  triangular 
square  matrix  B such  that 


BB 


T 


= A. 


(Lower  triangular  means  all  entries  above  the  main  diagonal  are  equal  to  0;  B^ 
means  the  transpose  of  B. ) A method  for  computing  such  a B will  be  given  below. 

Returning  to  the  main  discussion,  let  J2  denote  the  nxn  matrix 


It 


a -((®ik))- 


This  matrix  is  nonnegative  definite  because  it  is  the  correlation  matrix  for  the 
random  variables  tj(Pj,  t),  . . . , Tj(Pn,  t),  for  any  fixed  t.  Thus  there  exists  an 
nxn  matrix  B such  that 


T 

BB  — Q • 


Now  let  y j,  . . . , yQ  be  n independent  copies  of  a Gauss-Markov  process  with 
covariance  function  t— exp(-6  | tj ),  where  o2  is  the  appropriate  variance  and  1/6  is 
the  relaxation  time  discussed  in  Chapter  III.  Let  y denote  the  n-dimensional 
process  (y  j,  . . . , yQ)  and  let  be  defined  by 


?jT  = By  T. 


Then  the  desired  processes  t)  may  be  taken  as  the  n components  of  t}. 

That  is  because 


Eii<t)T  ij(s)  = E(By(t)T)(y(s)BT) 


= BE y (t)T  y (s)  BT 


= Bo2e-6|t"s|  IBT 


2 -6  t-8  _ 
= o e ' $2 


where  I is  .the  nxn  identity  matrix.  In  particular  each  has  covariance  function 
t-*o2  e~ 6 I4!  and  t|j(t)  and  Tik(t)  have  correlation  ojk  as  desired. 

In  practice  one  generally  needs  to  be  able  to  simulate  a discrete  version  of  the 

Gauss-Markov  process:  w(0),  w(A),  w(2A) for  some  suitably  small  value  of 

A > 0.  This  is  done  as  follows.  If  w should  have  covariance  function  t-*a2  e“^  1*1 
then  let  e (0),  e (A),  ....  be  a sequence  of  independent  draws  from  a Gaussian  dis- 
tribution with  mean  0 and  variance  o2.  Let  w be  defined  recursively  by: 


w(0)  = e (0), 
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and  for  k > 1, 


w(kA)  =e"<5Aw((k-l)A)  + /l-e‘2<5Ae  (kA). 


(V-l) 


Then  it  is  easy  to  check  that  this  w has  the  desired  properties. 

For  the  multivariate  process  rj  =(tj tj  ),  one  needs  e (kA)  to  be  n- 

dimensional  multivariate  Gaussian  with  covariance  matrix  a2I  and  one  takes 


*1  (kA)T  = e~6A  ((k-l)A )T  + /l-e_26A  Be  (kA)T. 


(V-2) 


Thus  we  are  now  in  a position  to  simulate  all  the  components  in  our  model 

for  stochastic  fluctuations  i,=lc  + rii  + Ti{,  i = 1 n.  The  q j and  | are  all 

independent  of  each  other  ana  of  the  r| and  the  r) . are  simulated  as  per  the  in- 
structions above.  No  explicit  recommendation  is  made  here  on  simulating  £ 
and  the  t] *;  see  reference  [a]  for  some  possibilities.  C 

For  moving  sensors  the  situation  is  only  slightly  more  complicated.  The 
average  velocity  V in  (III-l)  and  (III -2)  is  replaced  by  average  relative  velocity  so 
that  6 in  (V-2)  becomes  time  varying.  The  distance  rik  between  sensor  i and  k 
becomes  a time  varying  function  so  that  the  matrix  B also  becomes  time  varying. 

We  now  show  how,  given  a nonnegative  definite  $2,  one  can  compute  the  lower 
triangular  matrix  B such  that  BBT  = J2.  The  technique  is  called  the  Cholesky 
algorithm  (see  reference  [i],  p.  114). 

Suppose  = ((pik))  • Then  set  bjk=0  for  all  k>  i and  compute  the  other  bjk's 
by  the  following  algorithm.  (The  arrow  notation  means:  assign  the  value  at  the 
tail  of  the  arrow  to  the  variable  at  the  head  of  the  arrow. ) 

FOR  k = 1 TO  n 

, / k-1  2 

^kk*’  ^°kk  " jil  bkjy 

FOR  i = k + 1 TO  n 

bik  *"  ^ik‘  bkj^/bkk 

NEXT  i 
NEXT  k. 
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The  entire  procedure  can  be  summarized  as  follows: 

(i)  Given  the  propagation  loss  function  N^y  use  a Monte  Carlo  program 
to  estimate  the  autocorrelation  function  ^ over  a range  of  values 
which  includes  the  minimum  and  maximum  separation  of  the 
sensors.  Use  this  same  Monte  Carlo  to  estimate  the  mean  and 
variance  of  the  t|  process. 

(ii)  Compute  the  matrix  B such  that  BBT=  ((^(r^)))  where  r^  is 
the  distance  between  sensor  i and  sensor  k.  B is  used  as 
described  above  to  generate  the  processes  r^,  . . . , r)n. 

(iii)  For  a given  target  path,  compute  the  functions  x.(kA)  + £ .(kA) 
for  i = 1,  ...»  n,  and  for  k =0,  1,  . . . , until  some  inequality 
Xj(kA)  + i j(kA)  > 0 is  satisfied.  Store  this  value  kA  and  repeat 
the  procedure  building  an  empirical  distribution  function  which 
estimates  cdp. 

Note  that  the  computation  of  B is  a one-time  operation,  if  the  sensors  are 
stationary.  For  this  case  the  simulation  of  the  detection  process  is  only  trivially 
more  complex  than  it  is  if  we  assumed  that  all  sensors  were  mutually  independent: 
the  matrix  operation  B e must  be  performed  at  each  time  step  kA.  For  moving 
sensors  the  computations  in  step  (ii)  above  must  be  performed  at  every  time  step 
and  this  may  be  prohibitive  except  for  large  scale  computer  simulations. 


Symbol 

x 

I 

FOM 

NW 

cdp 


*c 

^i 

’•i 

pik 

Sj(t) 

j 


o 

j 

V. 

J 

Tl(P,t) 


d (P.Q) 

t(P.t) 

X 
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LIST  OF  SYMBOLS  USED 


Definition  Page  First  Used 

deterministic  component  of  signal  excess  1 

random  component  of  signal  excess  1 

sonar  figure  of  merit  1 

propagation  loss  function  1 

cumulative  detection  probability  1 

deterministic  component  of  signal  excess  2 

at  sensor  i 

random  component  of  signal  excess  at  sensor  i 2 

residual  component  of  % j,  . . . , £n  3 

correlated  component  of  % . 3 

independent  component  of  % . 3 

correlation  between  ^(t)  and  r^(t)  3 

location  of  j**1  noise  source  at  time  t 5 

index  set  for  the  noise  sources  5 

radiated  noise  from  the  j**1  noise  source  5 

initial  course  of  j1*1  noise  source  5 

speed  of  noise  source  5 

ambient  noise  at  point  P at  time  t 5 

Euclidean  distance  between  points  P and  Q 5 

the  power  process,  t)  = 10  log  £ 6 

intensity  of  Poisson  process  generating  the  initial  6 
positions  of  the  noise  sources 

finite  number  of  noise  sources  for  the  bounded  noise  6 
process 

region  in  which  noise  sources  are  distributed  for  the  6 
bounded  noise  prooess 
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LIST  OF  SYMBOLS  USED 


Symbol 


N 


<P 

Z, 


a 

10/3 

Tj(t) 

* 


v 

6 


iyp) 

^(u) 


rij 

n 

B 


Definition 


probability  distribution  for  the  {Vj} 
probability  distribution  for  the  {Nj} 
characteristic  function  of  £(P,  t) 

Z.  - loV10 

j 

f(r)  = 10  NW<r)/10,  power  loss  factor 
spreading  constant 


attenuation  coefficient 


distance  of  S^ (t)  from  fixed  point  P 
autocorrelation  function  of  the  time  series 

t “*  T)  (t) 

mean  of  the  speed  distribution,  / vdFy(v) 

exp(-6t)  is  approximately  the  autocorrelation 
function  for  the  time  series  t-»ri(t);  1/6  is 
called  the  relaxation  time 


distance  of  S.(t)  from  point  P at  fixed  time  t 

correlation  in  ambient  noise  between  points 
separated  by  distance  u 

distance  between  sensor  i and  sensor  j 


matrix  of  correlations 


square  root  of  n , BBT  = fi 


((^<v)) 
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